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^ ■ Abstract 

_>: 

'k>( \ We study the single spin asymmetry (SSA) for the pion production with large transverse 

V^ ■ momentum pf in semi- inclusive deep inelastic scattering ep^ — >■ evrX. We derive the twist- 

3 cross section formula for SSA, focussing on the soft-gluon-pole contributions associated 
with the twist-3 distribution for the nucleon and with the twist-3 fragmentation function for 
the pion. We present a simple estimate of the asymmetries due to each twist-3 effect from 
nucleon and pion, respectively, by fixing the overall strength of the relevant nonperturbative 
quantities by the data on the SSA A^ in p^p -^ nX collision. 
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1 Introduction 

Single transverse spin asymmetries (SSAs) in the strong interaction have a long history since 
70s and 80s when the large asymmetries were observed in pp -^ A^X [1] and p^p — > nX [2] 
in the forward direction. These data triggered lots of theoretical activities to clarify the 
mechanism of the asymmetries [3]. Experimentally measurements of SSA at higher energies 
in pp collisions have been performed at RHIC [4, 5, 6], and SSAs in semi-inclusive deep 
inelastic scattering (SIDIS) have been also reported [7, 8]. The SSA is a so-called "naively 
T-odd" observable proportional to {px q) ■ S±, and time-reversal invariance in QCD implies 
that it only occurs from the interference between the amplitudes which have different phases. 
In the literature, two QCD mechanisms have been used for describing the observed large 
SSAs. One is based on the use of so-called "T-odd" distribution and fragmentation functions 
which have explicit intrinsic transverse momentum k_i_ of partons inside hadrons [9, 10]. 
This mechanism describes the SSA as a leading twist effect in the region of small transverse 
momentum pt of the produced hadron. Factorization formula with the use of k_i_ dependent 
parton distribution functions and fragmentation functions has been extended to SIDIS [11], 
applying the method used for e"'"e~ annihilation [12] and the Drell-Yan process [13], and the 
universality property of those A;^-dependent functions have been examined in great details 
[14, 15, 16, 17, 18]. Phenomenological applications of the "T-odd" functions have been also 
performed to interprete the existing data for SSAs [19]. 

The other mechanism describes the SSA as a twist-3 effect in the coUinear factorization, 
and is suited for describing SSA in the large pt region [20, 21, 22, 23, 24, 25, 26]. In 
this framework, SSA is connected to particular quark-gluon correlation functions on the 
lightcone. This method has been extensively applied to SSA in pp collisions, such as direct 
photon production p^p -^ 7X [21], pion production p^p — >■ ttX [22, 23, 26], and the hyperon 
polarization pp —^ 1\}X [24] etc. Although the above two mechanisms describe SSA in 
different kinematic regions, it has been known that a soft-gluon pole function appearing 
in the twist-3 mechanism is connected to a A;_l moment of a "T-odd" function [16]. More 
recently, the authors of ref. [27] showed that the two mechanisms give the identical SSA 
for the Drell-Yan process in the intermediate Pt region, unifying the two mechanisms. 

In this paper, we study the SSA for large-pT pion production in SIDIS, ep^ — ^ cttX. ^ 
The QCD factorization tells us that the cross section for this process is generically ex- 
pressed as convolution of three quantities, a distribution function for the transversely po- 
larized nucleon, a fragmentation function for the pion, and a partonic hard cross section. 
In the framework of the collinear factorization, the SSAs at the leading order in QCD per- 
turbation theory require participation of an additional, coherent gluon into the partonic 
subprocess [20] , and the corresponding gluon field can be generated from either the nucleon 
or the pion, as shown in (a) and (b) of Fig. 1, respectively. These two contributions give 
rise to the following twist-3 cross section, denoting respectively as the (A) and (B) terms: 

(A) GF{xi,X2)®D{z)®dK, 
^The main result for SSA was reported in [28] by two of the present authors (H.E. and Y.K.). 



(B) 5q{x)(g)EFizi,Z2)(S)aB, (1) 

where Gf{xi,X2) and Ef{zi,Z2) are, respectively, the twist-3 distribution function for the 
transversely polarized nucleon and the twist-3 fragmentation function for the pion; the 
definition of these twist-3 functions is given in (2) and (8) below. In (1) there also appear 
the familiar twist-2 functions: 6q{x) is the quark transversity distribution of the nucleon, 
and D{z) is the usual quark/gluon fragmentation function for the pion. cta,b are the partonic 
hard cross sections, and we present their explicit formula in QCD perturbation theory. We 
note that both distribution and fragmentation functions in the (A) term of (1) are chiral- 
even, while both of those functions in the (B) term are chiral-odd. Accordingly, we have 
only one diagram in Fig. 1(b) relevant to the (B) term. We also note that there exist purely 
gluonic twist-3 distributions for the transversely polarized nucleon which can contribute to 
the (A) term in (1) [29]. The study of this term is beyond the scope of the present work, and 
we shall restrict ourselves to the kinematic region where the contribution from the quark- 
gluon correlation is expected to be dominant compared with that from the three-gluon 
correlation. 

It has been already shown that the twist-3 functions Gp and Ep appearing in (1) can 
potentially be important sources for the large asymmetry Ajv observed at large xp for 
p^p — > 7iX [22, 26]. Therefore, our interest here is to study how the above effects (A) 
and (B) contribute to the SSA in SIDIS, if the overall strengths of the relevant twist-3 
nonperturbative quantities Gf{xi,X2), Ef{zi,Z2), as well as the transversity Sq{x), are 
determined so as to reproduce A^v in pp collisions [2, 4]. 

The remainder of this paper is organized as follows: In section 2, we present definition of 
the twist-3 distribution and fragmentation functions necessary in our analysis. In section 3, 
after a brief description of the kinematics for SIDIS, we present the polarized cross section 
formula for the two contributions in (1), focussing on the soft-gluon-pole contributions. 
Azimuthal asymmetries from each contribution are discussed in detail and their numerical 
estimate is also presented. A brief summary is given in section 5. Appendix provides an 
explicit relation between the different definitions for the twist-3 distribution functions. 

2 Twist-3 quark-gluon correlation functions 

We first discuss the definition and basic properties of the twist-3 distributions for the 
transversely polarized nucleon, which are relevant to ep^ -^ enX. As shown in (a) of Fig. 1, 
these distributions actually represent the quark-gluon correlation inside the nucleon. For the 
transversely polarized nucleon with momentum p^ and spin Sj^, there are two independent 
twist-3 quark-gluon correlation functions which are defined as nucleon matrix element of 
nonlocal lightcone operator [21, 22, 30] 

dX f djj, 



= ^ Wn- e'^^-'-GF{xuX2) + ^^ Mh S1Gf{xu x^) + ■ ■ ■ , (2) 






(a) 



(b) 



Figure 1: (a) Contribution to ep^ — ^ evrX from the twist-3 distribution Gp] in the left (right) 
diagram, a quark (gluon) fragments into 7r+anything. (b) Contribution to ep^ -^ cttX 
from the twist-3 fragmentation function Ep- The white circles represent partonic hard 
scatterings. Mirror diagrams obtained by the interchange between left and right of the cut 
also contribute. 



where V' is the quark field, F"^ is the gluon field strength tensor, n'^ is the light-like vector 
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0) with p ■ n = 1, e 



apnS± 



e '^'^Pxn^S±u with eoi23 = 1, and the spin vector satisfies 



Sj_ = —1, S± ■ p = S± ■ n = 0. [fin, An] = P exp [ig /_(* dtn ■ A{tn)] represents the gauge-link 
which guarantees gauge invariance of the nonlocal lightcone operator, and "■ ■ ■" denotes 
twist-4 or higher-twist distributions. The correlation functions Gp{xi,X2), Gp{xi,X2) are 
defined as dimensionless, and the nucleon mass M^ in the RHS of (2) represents a natural 
scale for chiral-symmetry breaking. Each correlation function depends on the two variables 
Xi and X2, where xi and X2 — xi are the fractions of the lightcone momentum carried by 
the quark and the gluon, respectively, which are outgoing from the nucleon. From P- and 
T-invariance in QCD, one can show the symmetry properties 



Gp{xi,X2) = Gp{x2,Xi] 



Gp{xi,X2) = -Gp{x2,Xi), 



(3) 



and from hermiticity both Gp{xi,X2) and Gp{xi,X2) are real functions. These two correla- 
tion functions Gp{xi, X2) and Gp{xi, X2) constitute a complete set of the twist-3 quark-gluon 
correlation functions for the transversely polarized nucleon (see e.g. ref. [30] for a concise 
discussion) . 

In the literature, however, another set of twist-3 quark-gluon correlation functions are 
also used [21, 22]: Replacing the field strength gF°'^{fin)nj3 in the LHS in (2) by the trans- 
verse components of the covariant derivative, D'^{fin) = 9" — igA']_{nn), the same Lorentz 
decomposition in the RHS defines the two twist-3 distributions Gd{xi,X2) and Gd{xi,X2), 
in place of Gp{xi,X2) and Gp{xi,X2), respectively. Physical meaning as well as the basic 
properties of Gd{xi,X2) and Gd{xi,X2) are similar to Gp{xi,X2) and Gp{xi,X2), except 
that P- and T-invariance implies Gp,{xi,X2) = —G£){x2,Xi), Gd{xi,X2) = Ge,{x2,xi). 

Overcomplete set of four twist-3 correlation functions Gp{xi,X2), Gp{xi, X2), G c,{xi, X2) , 
and Gd{xi.,X2) for transversely polarized nucleon has been conveniently used on a case- 
by-case basis to express the cross section formula in the literature. Here we present an 



explicit relation between the "F-type" distributions {Gf{xi,X2),Gf{xi,X2)) and the "D- 
type" functions {Gd{xi,X2), Gd{xi,X2)). For this purpose, we refer to the operator identity 
with manifest gauge invariance, 

D°'{lj,n)[ij,n, \n] =i dt[^n, tn]gF"^{tn)ni3[tn, An] + [nn, \n]D°'{\n), (4) 



which relates (2) with the corresponding matrix element for the D-tjpe distributions. One 
then obtains 

Gd{Xi,X2) = P Gf{Xi,X2), (5) 

Xi — X2 

Gd{xi,X2) = 5{xi - X2)g{xi) + P Gf{xi,X2)-, (6) 

Xi - X2 

where P denotes the principal value, and 

^(^) = Yr I ^^''""(P S±mm, Xnh4S^ ■ D{Xn)^{Xn)\p S^) 

Mj\f J Ztx 

1 /■ rl\ 1"°° - 

+^7^ / ^^^'^ / '^^'<^' - ^) ^P ^^IV^(0)[0' An]75^(?F^-"(/in)[/in, An]^(An)b 5±), 

ZMl^ J 111 J~OD 

(7) 

withe(/i — A) = 29{fi — X) — l. In the RHS of (5), the term proportional to 5(a;i — 0:2) vanishes 
due to the anti-symmetry of G£,{xi, X2) as stated above. In Appendix, we demonstrate that 
g{x) can be reexpressed in terms of the F-type distributions Gf{xi,X2) and Gf{xi, X2), and 
the twist-2 quark helicity distribution Aq{x) of the nucleon. The relations (5) and (6) show 
that the Z^-type functions are more singular than the F-type functions at the soft gluon 
point Xi = X2, while they are proportional to each other for Xi 7^ X2. In our analysis 
for ep^ -^ enX below, we assume that there is no singularity in the F-type functions, in 
particular, that Gf{xi,X2) is finite at xi = X2 {Gf{xi,xi) = 0, due to (3)). 

Similarly to the distributions, one can also construct the twist-3 fragmentation function 
for the pion with momentum i as [25, 16, 26] 

x{7r{i)X\'ijjj{Xw)[Xw,fiw]gF'^^{fiw)wi3[fiw, ~oow]\0) 



^(75/7.).,-^'^"'"'^ 



where we introduced the light-like vector w by the relation i-w = 1, and "■ ■ ■" denotes twist- 
4 or higher, and we again use the nucleon mass Mj^ as a generic QCD mass scale in order to 



define Ef{zi,Z2) as dimensionless. A fragmentation function obtained by shifting the field 
strength gF"^^{jj,n)wp from the matrix element (7r(£)X| ■ ■ ■ |0) to (0| ■ ■ ■ |7r(£)X) is expressed 
in terms of Ep{z2,Zi). Unlike (3) for the twist-3 distributions, time reversal invariance 
does not bring any constraint on the symmetry properties of the twist-3 fragmentation 
functions. As in the case for the twist-3 distribution, we can also define another twist- 
3 quark-gluon fragmentation function £'^(-21, -22), replacing (yfF"^(/iw)w/3 by the covariant 
derivative D ^{fiw) in (8) [31]. The relation between Ed{zi, Z2) and Ef{zi, Z2) can be easily 
obtained by using the identity (4) as 

Ed{zi,Z2) =5[ ]e{zi) + P-—— -——Ef{zi,Z2), (9) 

v^i Z2/ U/^iJ ~ U/^2j 



where 



z 1 .^-^ f d\ 



e{z) = ^j^Y.j ^e-^/^(0|[-oot.,0]7>75^(0)|7r(£)X) 
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Here we emphasize again that Ed{z\, Z2) is more singular compared to Ef{zi, Z2) at zi = Z2, 
while they are proportional to each other for Zi 7^ Z2. As will be discussed in the next 
section, the cross section for ep^ -^ evrX receives important contribution from the "soft 
gluon point", Xi = X2 or Zi = Z2, of the twist-3 functions defined above. We thus use the 
corresponding F-type functions to express those contributions in the cross section. Other 
terms which receive contributions with xi ^ X2-, zi ^ Z2 are expressed in terms of either 
Z^-type or F-type functions without any subtlety. 

3 Single spin asymmetry for ep^ — > cttX 

3.1 Kinematics 

Here we present a brief description of the kinematics for the SIDIS, e{k) + p{p, S±) -^ 
e{k') + 7r(£) + X. (See refs. [32, 33] for the detail.) We have five independent Lorentz 
invariants, Sep,Xbj,Q'^,Zf, and q^. The center of mass energy squared, Sep, for the initial 
electron and the proton is 

Sep = ip + ky-2p-k, (11) 

ignoring masses. The conventional DIS variables are defined in terms of the virtual photon 
momentum q = k — k' as 

X,, = -^ , Q' = -q' = -{k - k'f. (12) 

2p ■ q 



For the final-state pion, we introduce the scahng variable 

% = ^ . (13) 

p- q 

Finally, we define the "transverse" component of q, which is orthogonal to both p and i: 

gf = q^- ^-^v^ - P^F . (14) 

p ■ t p ■ i 

qt is a space-like vector, and we denote its magnitude by 



Qt = \J-qt ■ (15) 

To completely specify the kinematics, we need to choose a reference frame. We shall work 
in the so-called hadron frame [32], which is the Breit frame of the virtual photon and the 
initial proton: 

q^^= (0,0,0,-g), (16) 

Further, in this frame the outgoing pion is taken to be in the xz plane: 

As one can see, the transverse momentum of the pion is in this frame given by ir = Zfqx- 
This is true for any frame in which the 3-momenta of the virtual photon and the initial 
proton are collinear. By introducing the angle between the hadron plane and the lepton 
plane, the lepton momentum can be parameterized as 

/c^ = — (cosh-?/', sinh-?/'cos0,sinh-?/'sin0, — 1) , (19) 



and one finds 






cosh^= ^,^ -1. (20) 

We parameterize the transverse spin vector of the initial proton S^ as 



5^ = (0, cos $5, sin $5,0), (21) 

where $5 represents the azimuthal angle of S± measured from the hadron plane. With the 
above definition, the cross section for ep^ — > enX can be expressed in terms of Sep, Xbj, Q^, 
-2/5 Qt^ 'P ^"^^ ^s in the hadron frame. Note that and $5 are invariant under boosts in 
the g^direction, so that the cross section presented below is the same in any frame where q 
and p are collinear. 



3.2 Twist-3 polarized cross section 

We now proceed to derive the cross section for ep^ -^ enX, applying the method developed 
in [21, 22]. We use Feynman gauge in the actual calculation. In this framework, the 
phase necessary for SSAs is provided from the partonic hard cross section. To illustrate, 
consider the diagrams in the first line of Fig. 2(a), which are some of the leading order 
contributions in perturbation theory corresponding to Fig. 1(a). These diagrams involve 
the coherent gluon from the nucleon, participating in the partonic hard scattering. After 
the collinear expansion in the transverse momenta kj_ of the partons inside nucleon to 
0{kj_), the "external" partons for the partonic hard cross section are collinear to their 
parent hadron. Denoting the momenta of the external quark lines on the LHS of the cut 
as xip and i/z, and the momentum of the coherent gluon as {x2 — xi)p, where xi, X2, and 
z are the corresponding momentum fractions to be integrated over, the parton propagator 
coupled to the coherent gluon gives a factor l/{xi — X2 + is). Then the integration for 
the parton momentum fraction produces the phase from the pole contribution at the soft 
gluon point Xi = X2 (soft-gluon-pole (SGP)) [21], which may be evaluated through the 
distribution identity 

1 1 

P i7T6{xi-X2). (22) 



a;i — a;2 + ie xi — X2 

It is straightforward to see that another propagator in the LHS of the cut in the first and 
third diagrams in the first line in Fig. 2(a) is proportional to l/{xi~ie) = P {1 / xi) +iTT6{xi) , 
which produces the soft-fermion-pole (SFP) contribution at xi = [21]. Similarly another 
propagator in the LHS of the cut in the second and fourth diagrams in the first line in 
Fig. 2(a) is proportional to l/{xi—Xbj+ie) = P{l/{xi—Xbj)}—in6{xi—Xbj), which produces 
the hard-pole (HP) contribution at xi = Xbj [34]. There are other leading order diagrams of 
the type of Fig. 1(a), which receive the SFP and HP contributions. As anticipated, the SGP 
contributions, as well as the SFP and HP contributions, are eventually associated with the 
quark-gluon correlation beyond twist-2, and are expressed by a complete set of the twist-3 
quark-gluon correlation functions, Gp{xi,X2) and Gf{xi,X2) of (2). From the symmetry 
(3) under xi ^^ X2, Gf{xi,X2) does not contribute to the SGP at Xi = X2- On the other 
hand, as a result of the collinear expansion, some of the SGP contributions appear with 
the derivative like ■^Gf{x,x), while SFP and HP contributions do not appear with such 
derivative.^ In the large x region, one has \x-^Gf{x,x)\ 3> \Gf{x,x)\, since Gf{x,x) is 
expected to behave as Gf{x,x) ~ (1 — x)^ (/9 > 0). Accordingly, it is justified to keep 
only the derivative terms with -^Gf{x, x) in the SGP contributions for the case where only 
large x region is relevant. 

Similarly, the diagrams in Fig. 2(b), which are some of the leading order contributions 
in perturbation theory corresponding to Fig. 1(b), also receive the pole contributions. After 

^The SFP and HP contributions can be straightforwardly calculated in the lightcone gauge A+ — 0, in 
which we can make direct translation A-^ — > iF^'^/{xi — X2), corresponding to the relations (5), (6) for 
xi ^ X2, and no collinear expansion is necessary. Thus no derivative terms appear for those contributions. 
This feature was also confirmed by explicit calculation in the Feynman gauge in a recent paper [43]. 





(b) 



Figure 2: (a) and (b) represent Feynman diagrams which give rise to SGP contributions 
for Gp and Ep-, respectively. ® denotes fragmentation insertions. In both cases, mirror 
diagrams also contribute. 



the collinear expansion, the quark momenta entering and leaving the pion fragmentation 
function can be parameterized as (!,/ zi and 1/ Z2-, respectively. In this case, the pole con- 
tributions, which give the phase for SSA, come from the SGP contribution at zi = Z2 and 
the HP contribution at Zi = Zf{l — Qt/Q'^)- But there appears no SFP contribution, since 
1/zi ^ {i = 1, 2). Again, some of the SGP contributions appear with the derivative of 
Ep, while no derivative appears for the HP contributions. ^ 

In our present study, we shall focus on the derivative terms of the SGP contributions 
following the previous studies [22, 23, 24, 26]. This approximation should be valid if either 
Xbj or Zf is large, where the main contribution is from large- a; and large- 2; regions (see the 
cross section formula derived below). Also, physically, one expects that the SFP and HP 
contributions are suppressed compared to the SGP contributions, because there may be 
lots of soft gluons in the hadrons; in particular, for the SFP contributions, not much soft 
fermions are expected to exist in a hadron. 

The diagrams which receive the SGP contributions of the type (A) in (1) are shown in 

"^In principle, the principal- value part of the distribution identity of the type (22), combined with the 
imaginary part of Ep, could produce SSA. Discussion on such effects associated with nonperturbative 
generation of the strong interaction phase is beyond the scope of this work. Correspondingly, in the 
following, we shall assume Ef{zi, Z2) is real. 



Fig. 2(a). By including only the derivative terms of the SGP contribution, one obtains the 
spin-dependent differential cross section as 

dxhjdQ'^dzfdq^dcj) Sirxl^S'^pQ^ \ Q"^ J ^^ ^ " Jxrm„ x A^™ z 

X (^GUx, x)\ [Da{z)Caa^^ + Dg{z)Cga^^) ^ ^ (5(x - xq), 

(23) 
where aem = e^/47r is the QED coupling constant, and we have introduced the variables 

x=^, z = ^, (24) 

X z 

Xo^x.(l + ^|). (26) 

In (23), Gp denotes the twist-3 distribution for the quark flavor a defined in (2), Da{z)/Dg{z) 
is the unpolarized quark (flavor a)/gluon fragmentation function for the pion, and Cq is the 
electric charge for the quark flavor a. The color factor Ca^g is defined as Cg = Cp — Nc/2 = 
— l/{2Nc) for quark and Cq = l/{2Nc) for anti-quark, and Cg = Nc/2 for gluon. The 
partonic hard cross section aj^ (j = q, g) in (23) is written as 

afq = [1 + cosh^ iIj\ ajq — 2a|g — sinh 2-?/' cos a|g + sinh^ ifj cos 20 a^^, (j = g, g) (27) 

where 






?1 0:C;;? J I ^ I //n2 ^2 



+ Q^-# 



^Jg = 2aJ^ = 8x2 



a,^, = 4x2--(g^ + g4), (2^ 



^T 



and 



a^ =2x(l-2)^^^ + il^I): g^-^^' 
^OT = 2^OT = 82(1 - z) , 
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We note that the above Bj^ {j = Q, 9] /c = 1, ■ ■ ■ , 4) for the twist-3 cross section (23) are the 
same as the corresponding unpolarized hard cross sections [35, 33]. The suppression factor 
{Mj^qx/Q^) in (23) characterizes twist-3 nature of the cross section. We also note that, in 
the present case where either x^j or Zf is large, the main contribution to (23) is from regions 
with large x and z, and this allows us to include only the contributions from the valence 
component of the distribution functions and the favored component for the fragmentation 
functions. 

Likewise we focus on the SGP contributions for the (B) term in (1), which come from 
the diagrams shown in Fig. 2(b). We keep only the corresponding derivative terms of Ep, 
following [22, 23, 24, 26]. The spin-dependent cross section is obtained to be 

dxbjdQ'^dzfdqj^dcp 8nxh S'^Q'^ 



f Mn7t\ V- 2 /"^ dx f'^ dz ^ , . 



X {2 



^^Im^u^) 



(1 - 2x)g2 / ^ d ^^, , 

xqj' + (1 — x)Q^ \ dz 



xCa^ 



X 



1 — X 



6{z- zq), 



(30) 



where Ep is the twist-3 fragmentation function for the quark flavor a defined in (8), the 
color factor Ca is defined as Cg = —1/{2N^ and Cq = l/{2Nc) for quark and anti-quark, 
respectively, the variable Zq is given by 



Zq = Zf\l + 



X qfj. 
and the partonic hard cross section a^ is defined as 



CT^ = (1 + cosh^ il)) sin $5 



-4ggr 



sinh2-?/'sin($5 — 



4g2 



£(g2 + q^) 



(31) 



+ sinh^ ip sin(<l>5 - 20) 



-4Q=^ 



SqriQ^ + Qt) 



(32) 



Here again we note the presence of the suppression factor M^/Q in (30) which characterizes 
the twist-3 cross section. Similarly to (23), with either x^j or Zf being large, we can restrict 
the sum over a to the contributions from the valence component of the distribution functions 
and the favored component for the fragmentation functions in (30). Note that, in (30), 



there are two terms with 



{d/dzi)EF{zi,z) 



not have definite symmetry property under zi 



and {d/dz)EF{z,z), since Ef{zi,Z2) does 

* Z2- 



3.3 Azimuthal asymmetry 

A remarkable feature of the polarized cross section for ep^ -^ enX is its characteristic 
azimuthal dependence. From (23) and (30), one finds the cross sections can be decomposed 



11 



as 






and 



dxbjdQ'^dzfdqj'C 



a^ sin($s) + fff sin($5 - 0) + o"^ sin(<l>5 - 20). (34) 



These should be compared with the similar decomposition of the twist-2 unpolarized cross 
section [35, 32]: 

^5_.unpol 

= a^ + a^ cos(0) + ff^ cos(20). (35) 



dxbjdQ'^dzfdq'j'dc 

As in (19) and (21), our azimuthal angles and $5 are defined in terms of the hadron plane. 
If one uses the lepton plane as a reference plane to define the azimuthal angle of the spin 
vector of the initial proton as 0^, and that of the hadron plane as 0^, as employed in [7], 
one has the relation $5 = 0/i — 05 and = 0/j. From this relation, one sees that azimuthal 
dependence of the a^ term in (33) is the same as the Sivers effect (oc sin(0/i — 0^)), and 
that of erf term in (34) is the same as the Collins effect (oc sin(0/j + 05)). Our cross sections 
(33), (34) have additional azimuthal components which are absent in the leading order cross 
section formula obtained with Sivers and Collins functions. 

We now proceed to estimate the SSA in SIDIS by using the obtained formulae (23) 
and (30). So far we don't have definite information on the nonperturbative functions, the 
twist-3 correlations Gp, Ep and the transversity distribution 6qa-^ In order to see the 
qualitative behavior of the SIDIS cross sections, we fix those nonperturbative functions 
so that the similar twist-3 mechanism associated with those functions can reproduce the 
observed An for p^p -^ nX [2, 4]; the twist-3 cross section for p^p -^ nX also receive 
the SGP contributions analogous to the (A) and (B) terms of (1), although the additional 
parton distribution associated with the unpolarized nucleon participates. For simplicity, we 
require that each of (A)- and (B)-type contributions can reproduce A^ for p^p -^ nX: The 
authors of ref. [22] used the ansatz Gp{x,x) = KaQaix) with a fiavor-dependent constant 
Ka to model the valence component of Gp{x,x) using the unpolarized quark distribution 
qa{x), and found that the SGP contributions associated with the derivative ■£:Gp{x,x) 
approximately reproduce An observed in [2, 4] with the value Ku = —Kd = 0.07. 

As was stated below (8), Ep{zi, Z2) does not have definite symmetry property. Here we 
assume 2[{d/dzi)Ep{zi, z)]^-^^^^ = {d/dz)Ep{z,z) for simphcity. With this assumption, the 
SGP contributions associated with the derivative oi Ep also reproduce An in [2, 4], using the 

^Ep{z,z) is related to the fc^-moment of the difference of the two types of the Cohins fragmentation 
functions with the "future-pointing" and "past-pointing" gauge hnks [16]. Ahhough it's been shown in 
[17] that one can derive factorization for SIDIS as weU as e'^e~ annihilation by using the "future-pointing" 
gauge link, no definite relation is known between the two types of the Collins functions with different gauge 
links. Therefore the above relation does not constrain Ep{z, z) at present. 
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d.nsd.tz Ef^iz^z) = KaDa{z) with K^ = K^ = -Kj = -Kl = -0.19 and combining 
with the assumption for the transversity, 5qd{x) = Aqd{x) and 6qu{x) = 0.58Ag„(x) [26]. ^ 
In our estimate, we use GRV unpolarized parton distribution qa{x) [36], GRSV polarized 
parton distribution Aqa{x) [37], and KKP pion fragmentation function Da,g{z) [38]. 

To see the magnitude of each component in (33) and (34), we calculate the azimuthal 
asymmetries normalized by the unpolarized cross section. We define ^-integrated azimuthal 
asymmetries as, setting $5 = 7r/2, 

(T (7 (7 

(l)w = -^. (cos0)iv = TTo, (cos20)jv = -%, (36) 

<7q 2a^ 2a^ 

for the (A) contribution, and 

R R R 

(1)^ ^ %, (sin(<|.5 - <P))n = ^, (sin($5 - 20))Ar = ^, (37) 

for the (B) contribution. For the purpose of illustration, we choose two sets of kinematic 
variables in the region where our formula is valid. The first one is Sep = 300 GeV^, 
Q^ = 100 GeV^, Xbj = 0.4, which is close to the COMPASS kinematics. Another one is 
Sep = 1000 GeV^, Q2 = iqo GeV^ and Xbj = 0.12, which is in the region of planned eRHIC 
experiment [39]. Both sets give the same cosh-?/) in (20). 

Fig. 3 shows the SSAs (36) for ep^ -^ ctt^X from the (A) contribution. One sees that 
(l)Ar can be as large as 5% at x^j = 0.4, while (cos0)7v and (cos 20) at are much less than 
1%. For comparison we also showed the azimuthal asymmetries for the unpolarized cross 
section (35), (cosn0)o = o"n/(2o"?) {n = 1, 2), in the same figure. (l)Ar is comparable 
to (cos0)o- Here the derivative {d/dx)Gp{x,x) together with Gp{x,x) = Kaqa{x), which 
explains rising Aj^ for p^p -^ nX at large xp, causes similar rising (l)Af toward Zf —^ 1. 
We also note the SSA has smooth g-p- dependence. (Recall that the transverse momentum 
of the pion, ip, is given by ip = Zfqx as (18).) 

Fig. 4 shows the azimuthal asymmetries (37) from the (B) contribution. At both energies 
Sep = 300 and 1000 GeV^, these asymmetries are much smaller than the (A) contribution 
of Fig. 3, although the derivative {d/dz)Ep{z,z) with Ep{z,z) = KaDa{z) causes rising 
behavior of SSA at large Zf. From these calculations, the effect of the (B) contribution is 
negligible in ep^ -^ evrX, even if one fixes the strength of the nonperturbative functions so 
that the corresponding (B)-type contribution solely reproduces A^ for p^p -^ nX. 

4 Summary 

In this paper, we have studied the SSA for the large-p^ pion production in SIDIS ep^ -^ 
cttX. We derived the twist-3 cross section formula, focussing on the the SGP contributions 

'"'In [26], we employed the same choice for the d-quark as above, but used K^ — —0.11 and 6qu{x) = 
Ag„(x) for the w-quark, instead of the above values. From the isospin and charge conjugation symmetry, 
however, one should have K^ = KJ — —K"^ — —K^ . Combined with the new choice for the w-quark 
transversity, Squ{x) — 0.58Ag„(x), the value K^ — —0.19 gives the same prediction for An as in [26]. 
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associated with the derivative of the twist-3 distribution Gp for the nucleon, and with the 
derivative of the twist-3 fragmentation function Ep for the pion. This approximation should 
be vahd in the region where xi,j or Zf is large. In estimating the impact of each contribution, 
we fixed the overall strengths of Gp and Ep so that each of them independently reproduces 
Atv in p^p — > ttX, and found that the Ep contribution is negligible compared to the Gp 
contribution. A full calculation including the nonderivative terms of the SGP contributions 
as well as the SFP and HP contributions will be reported elsewhere. 

Note added: 

After the submission of this work, the preprint [43] appeared. The authors of this 
reference have performed the twist-3 calculation for the gg-contribution (left figure of Fig. 
1(a)) to a^ term in (33), including all the SGP (both derivative and nonderivative terms) 
and HP contributions, and demonstrated its equivalence with the Sivers effect in the small 
Pt region. Their result for the derivative term of the SGP contribution agrees with ours. 
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Appendix: Relation between twist-3 distributions 

In the main text, we have discussed the relations (5)-(7) between the D-type and F-type 
twist-3 correlation functions for the transversely polarized nucleon. This appendix is de- 
voted to demonstrate that g{x) of (7) can be expressed in terms of the F-type twist-3 
correlation functions and the twist-2 quark helicity distibution Aq{x) for the nucleon. Such 
relation involves the components with different number of partons, and can be obtained as 
a consequense of constraints from the equations of motion and Lorentz invariance. 

We start with the well-known definition of the spin-dependent quark distributions of 
the nucleon [40, 30], 

^e^^^(p^|V5(0)[0, An]7'^75^(An)b^) = 2M^ [^q{x)p^{S ■ n) + gT{x)S^^] + ■ ■ ■ , (38) 

where Aq[x) and grix) are the twist-2 and -3 quark distributions for the longitudinally and 
transversely polarized nucleon, respectively, and "■ ■ ■" denotes Lorentz structure of twist 
higher than three. Using the QCD equations of motion, it can be shown [41] that grix) can 
be expressed in terms of the D-tjpe functions as 

gpix) = —— dx' Gd{x,x') + Gd{x' ix) + Go{x.,x') — Gi:,{x\x) . (39) 



14 



0.02 

0.00 
-0.02 
-0.04 
-0.06 
-0.08 

0.10 

0.05 

0.00 h 
-0.05 
-0.10 



<COS2(/)>o 



<COS(^>0 



<COS<^>N <C0S2(/)>N Sep=300 GeV 

Q=10 GeV 

Xbj=0.4 

Zf=0.2 




<COS</»o 




0.03 

0.02 

0.01 

0.00 

-0.01 

-0.02 

-0.03 

-0.04 



<C0S2(^>0 



<cos2(j»N 

<COS(/)>N 




<COS<j)>o 



(b) 



Sep=1000 GeV^ 
<1>N Q=10GeV 
Xbj=0.12 
Zf=0.2 



7 8 

qj (GeV) 



10 



0.02 

0.00 

-0.02 

-0.04 

-0.06 




<1>N 

Sep=1000 GeV^ 
Q=10 GeV 
qT=7 GeV 
Xbi=0.12 

_i_ I I 



(d) 



<C0S2(^>N 
<COS(^>N 



<COS(f»o 

I I L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Zf 



Figure 3: The azimuthal asymmetries for the Gp contribution to ep^ -^ en^X at two 
energies, (a) and (b) show the gr-dependence, while (c) and (d) show the Zf dependence. 



Substituting the relations (5) and (6) into (39) and using the symmetrty properties (3), we 
get 

"i ,Gf{x,x') + Gf{x,x') 



9t{x) = -—{g{x) + P j dx' 



, (40) 

X — X J 

We note that this result has been obtained restricting the quark and gluon fields strictly 
on the lightcone. As is well-known, however, similar result can be obtained using the 
operator product expansion, where the fields are not restricted on the lightcone and thus 
the constraints from Lorentz invariance are fully taken into account. Here it is convenient to 
empoy the nonlocal version of the operator product expansion using the operator dentity [42, 
30] 



d 



d 



^M I ^^-^(0)7%[0,^]V'(^) - ^V'(0)7^5[0,2;]^(z) 
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Figure 4: The azimuthal asymmetries for the Ep contribution to ep^ — > erc^X at two 
energies, (a) and (b) show the gy-dependence, while (c) and (d) show the Zf dependence. 



^' dtip{0)[0, tz] /|275 (t - i) gF"''{tz)z, - ^gF'^P{tz)z, \ [tz, z]^{z) 



+ 



^(0)75cr"'';2p[0, z]t P^Piz) - ij{0)ip-f,a'''zp[0, z]i;iz) 



(41) 



where z^ is not restricted on the hghtcone. This identity is exact up to the irrelevant terms, 
i.e., the terms proportional to quark mass, the twist-4 {0{z'^)) contributions, and the total 
derivatives. In the lightcone limit z^ — ^ 0, the nucleon matrix element of operators on both 
sides can be expressed in terms of the distribution functions defined in (38): 



{p 5|^(0)[0, z]Yl5^{z)\p S) = 2Mn I dxe-'^'P-' Ip"^ ^ 

-1 \ P - z 



^q{x) - grix) + 0{z^ 
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+ S'' grix) + 0{z^ 



+ 



(42) 



and the F-type distributions of (2); note that the matrix element of the last line of (41) van- 
ishes using the equations of motion. Thus the matrix element of (41) yields the differential 
equation 

dgrix) ... X 

X — ; h Aq[x) 

ax 



2 



dx' 

-1 X — x' 



'd_ 

dx 



d 



+ J^]Gf{x.x') + 



dx 



'd_ 
dx 



_d_ 
dx 



- Gf{x,x') 



(43) 



where we have used the symmetry relation (3). The solution of this equation with the 
boundary condition, grix) = for |a;| > 1, reads 



9t[x) 



/ -^ Aq{xi) + -P dx2 



Xi — X2 



dxi dx2 I 



d 



d 



dxi dx2 



Gp{xi, X2) 



(44) 



where — 1 < x < 1, and the first term on the RHS gives the Wandzura-Wilczek part for 
grix) [30]. It is straightforward to see that the moments J dxx"'gT{x) using this solution 
reproduce the results obtained by the standard local operator expansion. Comparing (40) 
with (44), we find that g{x) is expressed in terms of Gf{xi,X2), Gf{xi,X2), and Aq{x), 
and get 



g[x) = —X 



e{x) 



dxi 



2Aq{x^) 



Xl 



^ TD A {Gf{Xi,X2) , Gf{Xi,X2) 

PI dX2< h (3^1 - X2) -, TTT 

1 I Xl - X2 (Xl - X2Y 



x\ 



(45) 



Combined with (6), this result shows that ^^(3^1,3^2) contains the contribution from the 
twist-2 operators corresponding to the Wandzura-Wilczek part, and the genuine twist-3 
part of the D-type distributions Gd{xi,X2), Gd{,xi,X2) can be completely expressed in 
terms of the F-type distributions Gf{xi.,X2) and Gf{xi.,X2)- 
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